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We consider non centered vortices and their arrays in a cylindrically trapped Bose-Einstein 
condensate at zero temperature. We study the kinetic energy and the angular momentum per 
particle in the Thomas Fermi regime and their dependence on the distance of the vortices from the 
center of the trap. Using a perturbative approach with respect to the velocity-field of the vortices, we 
calculate to first order the frequency shift of the collective low-lying excitations due to the presence 
of an off-center vortex or a vortex array, and compare these results with predictions which would be 
obtained by the application of a simple sum-rule approach, previously found to be very successful 
for centered vortices. It turns out that the simple sum-rule approach fails for off-centered vortices. 
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I. INTRODUCTION 

Vortices in trapped Bose-condensed gases have recently been observed at JILA and at ENS At ENS they 

have been formed by stirring the condensate with a focused laser beam with an angular frequency £1. Experimentally 
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there exists a threshold of the angular frequency of the stirring beam (fi c ) to nucleate a single vortex. When ~ fi c 
one vortex is created at the stable position at the center of the trap. But at Q > fl c and depending on the frequency 
Q*\ ' of the stirring beam it is possible to create configurations with different number of vortices forming vortex arrays 
CNJ . which are stable compared to a single vortex with correspondingly larger circulation. 

The presence of vortices in trapped condensates has been revealed by time of flight analysis Jl],^| and recently, by 
exciting the quadrupole oscillations of a stirred Bose-Einstein condensate ^ and observing their frequency splitting. 
With this method it has been also possible to measure the angular momentum induced by the rotation of the trap. 
Direct measurement of the angular momentum is difficult in atomic gases, since most of the diagnostic techniques 
provide information on the density profiles. As expected, for stirring frequencies below fl c no angular momentum is 
observed since no vortex can be nucleated, and just above f2 c the angular momentum has a jump of ~ h which 
indicates the presence of a stable singly quantized vortex at the center of the trap. However, for higher stirring 
frequencies, an apriori unexpected behavior of the angular momentum has been experimentally found ||: for ft > il c 
when vortex arrays are nucleated, the measured angular momentum does not show any jump of size fi as a new singly 
quantized vortex is nucleated, as is e.g. observed for the corresponding experiment in superfluid Helium Q; instead 
it is a smooth and increasing function of Vt. We shall see that this is not in contradiction with the quantization of the 
circulation of a vortex, but it is a consequence of the nucleation of off-axis vortices that form a vortex array. 

The purpose of the present paper is to provide a detailed analysis of the dependence of the kinetic energy and the 
angular momentum on the distance of an off-centered vortex with respect to the symmetry axis. We will generalize 
the results to vortex arrays to give a physical insight of the smooth behavior of the angular momentum found 
experimentally. 

In an axially symmetric trap the axial component of angular momentum (m) is a good quantum number, and in the 
absence of vortices the collective excitations of the condensate are degenerate with respect to the sign of m. However, 
the presence of vortices breaks this degeneracy, and elementary excitations carrying opposite angular momentum are 
no longer degenerate. The frequency shifts of the quadrupole oscillations due to the presence of a vortex on the 
symmetry-axis of the trap have been calculated in [gj using a sum rule approach. Analytic results for the energy 
splitting are obtained in in the Thomas-Fermi limit for the quadrupole modes and have been used in || to measure 
the angular momentum of the condensate. General expressions based on the hydrodynamic approach for the energy 
splitting due to a centered vortex have been obtained in |||7| within perturbation theory. 

In this paper we consider off-centered quantized vortex lines in large cylindrically confined condensates at zero 
temperature, and therefore, no dissipation mechanism || is taken into account. We calculate the kinetic energy 
and the angular momentum with respect to the symmetry axis |9[[Tc[| of an array of vortex lines in a trapped Bose 
condensed gas. Using the perturbative approach proposed in [Q and assuming the Thomas-Fermi limit, we generalize 
the analytic expressions for the energy splitting of the low-lying collective modes to consider the effect of the distance 
of off-centered vortices with respect to the symmetry axis. We also compare with the results obtained by a direct 



1 



extension of the sum-rule approach to off-center vortices, but find that such an extension, at least in the direct version 
examined here, fails. 



II. VORTEX STATES 



We consider a weakly interacting Bose-condensed gas confined in a harmonic trap T4 x t(r) at zero temperature. The 
condensate wave function can be written in terms of its density p(r) = ^(r)! 2 and phase S(r) by 

^(r) = ^pJrjexp[iS{r)]. (1) 

The number of atoms in the condensate is J dr\^f\ 2 = N, and the superfluid velocity is given by v = (h/M)VS, where 
M is the atomic mass. The ground state wave function in the absence of vortices has a spatially constant phase and 
therefore zero velocity, but when a quantized vortex is present the phase and the velocity field of the ground state 
have to be determined according to the boundary conditions of the system. 



A. Condensate Density 



The ground state wave function \&(r) is a stationary solution of the Gross-Pitaevskii equation 

^ 2 V 2 



2M 



+ Vext(r) + g | *(r) | 2 tf (r) = /i*(r) , (2) 



where the coupling constant g is given by the s-wave scattering length a through g = 4nh 2 a/M, and p is the chemical 
potential fixed by the normalization condition of the ground state. 

Non- rotating experimental traps have axial symmetry, with different radial (u>±) and axial (u> z ) trapping 

frequencies, whose ratio defines the anisotropy parameter A = oj z /uj±. So far, vortex arrays have been produced at 
ENS |2||| in a highly anisotropic cigar-shaped trap (quasi-cylindrically symmetric) with A ~ 0.05. Thus, for simplicity, 
we consider an idealized cylindrical trap that is uniform in the z-direction (ui z — 0) and has a harmonic confining 
potential in the radial direction of the form 

V ext {r x ) = ~Mujlrl, (3) 

with r\ = x 2 +y 2 . The harmonic trap frequency luj_ provides a typical length scale for the system, a± — (H/Muj^) 1 ^ 2 . 

We consider the Thomas- Fermi (TF) regime |Tl| ], valid for large condensates with positive scattering length, where 
the kinetic pressure can be neglected compared to the interaction energy-density. The ground state density of the 
condensate, in the absence of vortices, is given by 

^-fC 1 -^)- (4) 

for r± < R±, and po(r±) = elsewhere. The Thomas-Fermi radius of the cylindrical condensate is R± = 
(2p/Muj 2 _) 1 ^ 2 . The density is defined in the interval — R z < z < R z , where 2R Z is the length of the cylinder. 
In the cylindrical geometry the validity of the Thomas-Fermi approximation for the ground state is guaranteed by 
the condition Na/(2R Z ) — N±a 1, or equivalently, /i 3> hu>±, where N± — N/(2R Z ) is the number of atoms per 
unit length 0. From Eq. (|) 

N± = [ R± po 2nr ± dr ± = . (5) 
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When a quantized vortex is present at the position r^o, the density of the system drops to zero at the center of the 
vortex core whose size is determined by the healing length £. In the limit of large systems it can be approximated 
by £ = (^npoir^a)^ 1 ^ 2 , where po(r±o) is the density of the condensate at the position of the vortex core but in the 
absence of vorticity. For a centered vortex, Eq. ([|) yields po(0) = p/g, and the corresponding healing length £o = £(0) 
can be rewritten as 
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(6) 



implying the following set of TF inequalities for the length scales £o <C a ± *C Analogously, from (Q) the local 
healing length [j^J is 

Within the TF approximation, analytical expressions for the density of a vortex state have been obtained for a centered 
vortex HQ and for a straight off-axis vortex line [^3| . However, in weighted spatial averages of quantities varying only 
on the scale R± , the density of a vortex state may be replaced by the density of the vortex- free state, the corrections 
being only of order (£/R±) 2 = (a±/R±) 4 0, which is negligible. 



B. Velocity Field 

We consider states having a quantized vortex line along the z-axis and all the atoms flowing around it with 
quantized circulation. The property of a single-valued wave function leads to the quantization of the circulation 
around an arbitrary closed loop which encloses the vortex core 
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where the integer number k is the quantum of circulation. For sufficiently large frequencies Q > fl c a single vortex 
line with an integer quantum number k > 1 can, in principle, appear but this state is unstable and fragments into 
a vortex array formed by k vortices each with a unit of circulation and position r± Qi (1 < i < k) relative to the 
symmetry axis. 

We will not consider large arrays of vortex lines, where an average vortex density can be defined but small 

ones, corresponding to the experimental region [^]|| of multiple vortices (n < 5). Large vortex arrays may enter the 
" turbulent" region found experimentally. For small vortex arrays the density of the system can be approximated by 
the density of a vortex-free condensate, because vortices rotating in the same direction experience an effective repulsive 
interaction Q , and the inhomogeneities due to the well separated vortex cores are negligible in the condensate density 

Let us assume a generic quantized vortex line parallel to the z-axis at the position r^ an d with quantized circulation 
K. Then, from Eq. (0), Vxv = (h,K/M)6(r± — r± )z, where z is the unit vector in the z-direction, and the superfluid 
is irrotational everywhere except for the vortex core at rj_ = r±Q, where the density vanishes. 

For a positively oriented vortex in an infinite, uniform system, the velocity then is 



In particular, for a centered vortex only the tangential component carries non vanishing atomic flux, and v = 
Tik/ (Mr±)(p, where <p is the unit vector in the tangential direction in cylindrical coordinates (r±, cp, z). Eq. (||) once 
again shows that a cut-off of the order of the local healing length £ for the distance |r^ — r^ol is needed, because at 
distances smaller than y/2£ the velocity ([)]) surpasses the local velocity of sound c s (r±) = y / g\^j(r±)\ 2 /AI. 

For a confined system, the velocity field is affected by the boundary of the system and by the spatially varying 
density, and has to fulfill the following physical conditions: i) the normal velocity has to vanish at the boundary, 
and ii) the condition for stationary flow V • (pv) = has to be satisfied. It is well known from the rotating bucket 
experiment in superfluid helium |j,[l^,|l2| that for an homogeneous system confined in a cylinder of radius R±, the 
normal velocity vanishes at the boundary by introducing an oppositely oriented image vortex at r^i = (i?^/r^o) 2 r^o- 
The resulting velocity field is 



kzx(n-rio) kzx(n-rn) 



If the system has a density gradient, the condition ii) for stationary flow is fulfilled by introducing a small correction 
to the velocity which can be neglected when the density varies over a larger scale than the healing length |l2). Thus, 
in the TF limit the velocity field can be approximated as v ~ vo . 
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A vortex at r^o is influenced by the velocity-field induced by its mirror-vortex. If rj_o approaches the boundary 
this velocity diverges, signalling the break-down of the Thomas-Fermi approximation. A suitable lower cut-off for 
the distance from the boundary, ^boundary = -R±(£o/2-R_i_) 2 / 3 , is implied by restricting the velocities induced by the 
mirror-vortices by the local velocity of sound. 

From Eq. (10) one can see that the contribution of a vortex (vortex and image vortex) is additive. Therefore, the 
velocity field corresponding to a vortex array formed by n singly quantized vortices at rj_oi with i — 1...K is then 

K 

v(rj = ]>>(r ± ) (11) 

with 

/ h z x (t_l - rxoi) _ h z x (rj_ - r ± u ) 
\M \r ± -r ±Ql \ 2 M \r± - r X i 



where r±u = (R±/r±oi) 2 r±oi is the position of the vortex image corresponding to the i- vortex. 



III. KINETIC ENERGY 



In the Thomas-Fermi limit the dominant part of the excess energy of a vortex state over the ground state without 
a vortex is given by the kinetic energy of the velocity field H . It is obtained from (O.O) in the form 



Ekin = 2_j ^ i (^) 
i—1 i—1 j — 1 



with 



M f Rz f R± 
Ei = — dz dr±r±p (r±, z) 
1 J -r, Jo 
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d<pv 2 (r ± , Zl <p) (14) 



Eij = M I dz dr±r±p (r±,z) 
-Rz JO 

d(pvi(r±,z,ip)-Vj(rx,z,ip), (15) 

o 

where Ei is the self energy of the i- vortex and J3y the binary interaction energy between a couple of vortices i and 

j. Let us define the dimcnsionless quantities Ei — Ei/ ^ 2 ^ , and analogously E^ and Ekm- For an array of 

vortex lines along the z-axis in the idealized cylindrical trap described by Eqs. (||J), the integrals in ( [l4] , [l5| ) can be 
performed and we obtain after some calculations 1 , defining the rescaled distance- vector of the vortex-line from the 
z-axis x< = r± 0i /Rj_, 

E i = M 1 - *i) + (1 - *i) ln (^") (16) 

1 




2 2 9 

j ^i^j 



2x\x) 



x It is useful for this purpose to adopt a complex representation of the transverse vector rx — > z — x + iy and of the velocity 
field of a vortex at zo in the form Vj — > Vj = (h/M)iz/(r 2 L — zzq) with complex conjugate v* = —{%/M)i/(z — Zo) and to 
perform the integral over ip at fixed value of r±_ — \z\ as a complex contour-integration making use of the residuum-theorem. 
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(17) 



t j 

arctan [ -!— ^ iL ] . 

V 1 - X 4 • Xj ) 



Here the healing length \/2£o / is used as a cut-off. It is of interest to remark that for R± and £o fixed, both Ei and 
Eij are proportional to N . The expression ( [Tg| ) for the kinetic energy of a vortex in the presence of its image-vortex 
(self energy) agrees with the result presented by Fetter |14[ fo r one vortex. The kinetic energy of an array of vortices 
without their images was calculated by Castin and Dum [|15|. Due to the presence of image vortices, introduced to 
ensure the vanishing of the normal component of the velocity at the boundary, our result differs from that in (Tq | but 
it is only slightly more complicated. 

In Figure 1 we plot the kinetic energy of symmetrical arrays of singly quantized 1, 2, 3, 4 vortices, analogous to 
the experimental configurations of Ref. ||, as a function of the normalized distance from the trap axis. The kinetic 
energy goes to zero as the distance of the vortices from the boundary of the condensate goes to zero. This is a 
consequence of the boundary condition on the surface, where each vortex meets its image and thus is annihilated in 
the process. As mentioned in the preceding section a lower cut-off must be imposed on the distance of the vortex-line 
from the boundary. This is tantamount to restricting Xi by Xi < 1 — (Co/2i?±) 2 ^ 3 . The energy is maximal at the 
center of the trap, where the expression ( |l7| ) diverges logarithmically unless a cut-off on the minimal distance of the 
vortex-cores from the z-axis is imposed, whose size is again taken as V2£o/R±. Fig. 1 makes it clear that a single 
vortex or a vortex array must be created at the boundary where the energy required to set up the necessary velocity 
field becomes vanishingly small. However, as we have already mentioned, the description of a vortex close to the 
boundary is outside the scope of the present approach, because there the Thomas-Fermi approximation, on which it 
is based, breaks down. 



IV. ANGULAR MOMENTUM 



Let us recall the expression of the angular momentum of a vortex state in an imperfect Bose gas |9|,|l0[] . We consider 
an axially symmetric condensate and are interested in the angular momentum around its symmetry axis defined along 
the z-axis. For a single vortex at a distance r±o from the symmetry axis, the angular momentum of the system with 
respect to the z-axis must be calculated by integrating over the whole volume of the system: 



< L 



>= J d 3 rp{r) (r x v) , (18) 



where v is the superfluid velocity around the vortex and p(r) is the condensate density. As discussed already, we can 
neglect in the TF limit the effect of the vortex on the density profile and take p(r) ~ po( r )- Then, using cylindrical 
coordinates around the z-axis, we have po( r ) = po(r±, z )- I n the present geometry, L is parallel to the axis of rotation 
and its magnitude is given by 

<L Z >= J d 3 rp (r ± , z)z • (r x v) . (19) 

The triple scalar product can be written as z • (r x v) = r j_ v v , where v v is the tangential component of the velocity. 
Assuming that the boundaries of the system are given by — R z < z < R z , < rj_ < R±, then 

dz I p (r ± ,z)rj_dr 1 _ / v v ,r±dp. (20) 

-Rz JO JO 

It can be easily seen that, for fixed r±_ and z, the angular integral is a line integral around a closed path, which 



corresponds to the quantization of the circulation (|8j). Therefore, the angular integral in (20) will contribute only 
when the closed contour encloses the vortex core. That is, 



I 
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v v r±dip = (h v • dl = — 2ir 9(r± — r± a ) , (21) 
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where 8(r± — r±o) is the step function. Remarkably, the detailed form of the velocity field around the vortex drops 
out in this expression. In particular, the mirror vortex contributing in ([l2]) docs not contribute to this integral at all 
because it is positioned outside the condensate and therefore never enclosed by the integration contour in (£]]). This 
leads to the following result for the angular momentum 



< L, >= 



R, rR± 

dz I 



r_L Po(r±,z)dr A 



(22) 



Due to the step function, only the part of the Bose condensate at r± > r± , that is outside the smallest circle around 
the z-axis which still encloses the vortex-line, contributes to the total angular momentum and therefore to the moment 
of inertia of the condensate.^ 

We will consider a large cylindrical condensate, which provides a good approximation for very elongated cigar- 
shaped traps in the Thomas-Fermi regime. Inserting in Eq. ( P4 ) the Thomas-Fermi density (Q) it follows that the 
angular momentum per particle is 



< l z >= 



<L Z > 
N 



Tin 



R± 



-1 2 



(23) 



For fixed R± it is independent of TV. It is worth stressing that Eq. ( p3| ) is the dominant term of the angular momentum, 
where corrections of order (£/i?^) 2 = (a±/R±) 4 are neglected. Eq. ( p3| ) shows that even though the circulation (g) 
around each vortex is quantized in units of Ti/M the angular momentum per particle around the z-axis is not quantized 
in units of h, in general. 

Since L is linear in the velocity field, contributions from additional vortices are additive and the angular momentum 
per particle can be generalized to a vortex array of k singly quantized vortices at rj_oi {i — !...«) as 



< l z >= 



R± 



(24) 



It is interesting to note that for a singly quantized vortex (k = 1), < l z > is equal to Ti only when the vortex line is 
centered at the z-axis. Otherwise < l z > is strictly less than h and decreases when the position of the vortex core 
moves away from the center and approaches the boundary of the condensate. When the vortex core reaches the edge 
of the condensate (r±o — R±) then < l z >~ 0. Analogously, when k singly quantized vortices are present, the angular 
momentum per particle is lower than kH unless all cores are along the z-axis || . But this configuration corresponds 
to one vortex with circulation kTl/M which is unstable when n > 1 p7|| and breaks into an array of k vortices all with 
unit quantization fi/M. 

Let us now calculate the angular momentum per particle in the simplest vortex-arrays in cylindrical traps rotating 
with angular velocity f2. They are created if the array is permitted to reach a state of relative equilibrium minimizing 
the energy 



E w «(fi) =E kin (0)-Q < L z > 



(25) 



Here -Efcm(O) is the kinetic energy ( |13[ ) of the vortex-array in the non-rotating trap. As before we consider symmetrical 
arrays of 1-4 vortices. In Figure 2 the kinetic energy (|2^) is plotted against the common normalized distance of the 
vortices from the z-axis, for an angular velocity Q/2-k — 40 Hz, which is, in our idealized 2-dimcnsional trap and 
within the Thomas-Fermi approximation, the critical frequency where the 2-vortex array first becomes stable. The 
corresponding critical frequencies for the symmetrical 1, 3, and 4 vortex arrays are, respectively f2 c = 33.6,43.3, and 
46 Hz, where we have chosen for the sake of concreteness, the values of M, R± according to the experiment described 



2 The expression ( P4 ) can be generalized to take into account a more realistic three dimensional vortex line. Real axially 
symmetric traps are not cylindrical, but provide also a longitudinal harmonic confinement 7^ 0) . It induces an inhomogeneity 
and a density gradient of the condensate along the longitudinal direction which will affect the velocity. The resulting vortex 
line must then deform along its length |l5l|l6| and meet the boundary of the condensate at a right angle in order to satisfy 
the physical conditions for the velocity field. Therefore, the distance vector of a quantized vortex from the z-axis in general 
depends on z: r±o = r±o(z). By taking as a limit of radial integration r±o(z) in (0), the angular momentum expression ( |2^ ) 
allows to take into account such more complicated vortex configurations. 
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in (a]. As expected, the critical frequency is larger for large vortex arrays. These values are different from and much 
lower than the measured critical rotation frequencies at which vortices or vortex arrays are first observed to appear, 
because of the existence of energy-barriers which must be overcome before the relative energy-minima formed by the 
vortex states can be reached. 

It can be seen that in the rotating trap, for sufficiently high values of Q, a single vortex is in relative or absolute 
equilibrium only at the rotation axis as is of course well-known, while the symmetrical arrays of vortices have equilibria 
at finite distances from the z-axis which increase with the number of vortices because of their mutual repulsion. In 
Figure 3 we plot for the same vortex-arrays at the corresponding equilibrium configuration the average angular- 
momentum per particle as a function of f2. With increasing Q the equilibrium-positions of the vortices in the arrays 
move towards the z-axis and the average angular-momentum therefore increases. For a single vortex, since the stable 
position is at the trap center independently of the rotation frequency, the angular momentum is always lh. 

Let us remark here that, in the absence of dissipation, vortices created experimentally in rotating traps, in general, 
need not correspond to minima of (^5|). A dissipative mechanism must be active on the scale of the time- interval 
during which the rotation-frequency is switched on for their positions to be able to relax to an energy minimum. After 
the rotation is switched off the same mechanism will tend to lead to a relaxation of rj_o towards the minimum at the 
boundary r±o = R± of the energy (|l3|), which increases the observed value of r±Q. In the experiments on vortices 
in traps reported so far dissipation seems to play a negligible role. In the absence of such a mechanism, however, 
the vortices with distances r±a outside the energy-minimum cannot relax but experience a force and a corresponding 
Magnus deflection already in the rotating trap which leads to a rotation of an array of equidistant vortices around the 
z-axis even in the frame rotating with the trap. Vortices which are created at the boundary and have not yet reached 
their equilibrium-distance from the z-axis will then be observed to have an angular momentum which is smaller than 
the equilibrium value shown in Fig. 3. From this point of view, the fact that the first vortex which is formed is 
always observed to be in the center seems to indicate that the first vortex is, in fact, not created by the motion of 
a vortex-line from the boundary to the center but by a different mechanism, like e.g. the condensation of collective 
excitations with I = 1 into a vortex state with 1 = 1. 

Let us now turn to a discussion of some related experimental results which have been obtained in Refs. [§U^]. From 
the transverse absorption images of Ref. Q we estimate the angular momentum corresponding to the experimental 
configurations of the condensate with 1 up to 4 vortices. We proceed as in ||, first we obtain a qualitative measure of 
the ratio between the distance of each vortex from the center and the average radius of the expanding condensate, by 
measuring them in the portrayed images of the condensate after the time of flight. During the expansion the transverse 
lengths scale by the same factor , therefore the relation between the distance of a vortex from the center and the 
average radius of the condensate (r±o/R±) will be the same as before the expansion. From the transverse absorption 
images it can also be seen, that vortex arrays have uniform spatial distribution and thus, vortices are equidistant 
from the center as has been assumed in Figs. 2 and 3. Then, with the values of r±o/R± extracted from Q by using 
Eq. ( p^ ) the following estimates for the angular momentum per particle are obtained: < l z > /% = 1, 1.33, 1.36 and 
1.38, corresponding to the experimental configurations with 1, 2, 3 and 4 vortices, respectively. These values of 
< l z > /h are smaller than the equilibrium values given in Fig. 3 which seems to indicate that the relative radii of the 
vortex-arrays measured in j^] are larger than their calculated equilibrium values in the rotating trap. It would be nice 
to compare in detail the values of < l z > /fi we calculate from the images in [pi with those measured directly in jj). 
Unfortunately, however, the experimental parameters are slightly different in || and || which makes a comparison 
difficult. But it is worth stressing that the calculated result for the angular momentum qualitatively agrees with the 
experimental data in [§!: after the first jump of H corresponding to the nucleation of a singly quantized vortex centered 
in the trap, the angular momentum increases continuously when the number of vortices in the array increases [tHI , 
(which is equivalent to have increased the stirring frequency) without presenting other jumps of order Ti. 



V. ENERGY SPLITTING 



The presence of a quantized vortex in a confined condensate breaks time reversal symmetry: this produces a 
frequency shift of the modes with azimuthal quantum number ±|m|, which in the absence of vortices are degenerated. 
We want to calculate the frequency splitting of the low lying modes in large systems due to an off-center vortex. 

Theoretical calculations of the frequency shift produced by a centered vortex have been performed within different 
approaches: a sum rule approach J5j, a semiclassical approach based on a large N expansion Q, a hydrodynamic 
approach , and a full numerical solution of the linearized equations of motion |]l9| . 

Current experiments in ENS Q have excited the two transverse quadrupole modes m = ±2 of a quasi-cylindrically 
symmetric condensate. Due to the presence of a vortex, there exists a lift of degeneracy between the frequencies of 
these two quadrupole modes that causes a precession of the system. Measuring the frequency of precession and from 
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the analytical expressions of the frequency shift [H the angular momentum of the system is inferred. 

For simplicity, we will assume a large condensate trapped in a cylindrically symmetric trap ([|) and we will study 
the energy splitting for the low-lying modes due to an off-center vortex at a distance r$ from the symmetry axis. In 
the following we shall use the perturbative approach . (The perturbative result for the frequency-splitting of the 
to = ±2 modes induced by an off-centered vortex in a trap with different symmetry from the one considered here was 
reported in [p0[ ). Linearizing the Gross-Pitaevskii equation around the condensate using the decomposition ([!]) of W 
in condensate density and phase, and assuming the Thomas-Fermi approximation and the long-wavelength limit, it 
follows the coupled equations 

iu5p= V-(v Sp) + Av-OaoVSS) (26) 
M 

ioj5S=v Q -V5S+^5p (27) 

where 8p and SS are small deviations from the equilibrium values of the density and phase, respectively. Solving (|2 
for SS to first order in v and inserting in (p^) to eliminate SS, yields the perturbed wave-equation 



W^V-poV) *p = -*«V.(v <yp) + ^V-(poV(v -V*p)). (28) 



In first-order perturbation-theory the shift of the eigenvalue uj 2 is obtained from the expectation- value of the pertur- 
bation, taken with the unperturbed solution of (pSj), uj 2 , with vo = 0: 

co 2 -wg = -iuo < v ■ V > < (V-/9 V)(v -V) > (29) 

The expectation-value can be evaluated using the unperturbed Thomas-Fermi density (|j) if care is taken to interpret 
the radial part of the space-integral as a principal- value integral at rj_ = tq, because the actual condensate density at 
the position of the vortex strictly vanishes. 

Let us evaluate then the expectation-value for a single off-centered vortex ( |To| ) in the particular case of surface 
modes, for which fl2]J 



^=V^^+T e ™ V ' >*> = V*iu± (30) 

where we take to first as positive. We obtain after some calculation (and using to considerable advantage the 
technique referred to in the preceding footnote) the following result for the frequency shift, due to the off-centered 
vortex, exhibited by the surface mode with positive m 



(m+ 1) ' 
MR 2 V R 



-0 = ^^(1-^), (31) 



It turns out during the calculation that the mirror- vortex makes actually no contribution to the frequency-shift. 
Because the velocity- field enters linearly, the frequency shift induced by several vortices is additive. Replacing m by 
—777 is equivalent to inverting vo, i.e. uj — ujq simply changes sign. The frequency splitting between the two frequencies 
(with 777 and —to) is then 



1™ / 2m-2 

2m I _ r, 



»+-»- = mwA'-w^j- (32) 

It can be seen from (|3l| ) that there is no frequency shift for the dipole mode m = 1, at least not to first order in vo, 
i.e., the Kohn theorem, which predicts u)(m = 1) = lo±, is duely respected. 

The splitting of the quadrupole excitations (m = 2) has been measured in the ENS experiment j|. For this case it 
follows from (|32|) 



MR 2 



Ri 



(33) 



From Eq. (|33|) one can see that there is a term in the frequency shift that depends on the distance of the off-center 
vortex from the center. It means that the precession of the eigenaxis of the quadrupole mode, 9 = — ujj)/2\m\ 
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that is measured in ENS experiment |3j], will be slightly different for an off-center than for a centered vortex. And 
thus, from this difference the distance of the off-centered vortex can be inferred. 

It is worth noting that Eq. ([33]) in the particular case of a centered vortex (rn / R± = 0) does not give the same 
dependence on the trap and condensate parameters as the result obtained in because the symmetry of the trap 

considered in both calculations is different. 

The preceding calculation can be extended to axial vortex-modes |f22| travelling along the axis of the trap with 
wave-number k. For small k the unperturbed modes without vortex, to leading order in (kR±) 2 , take the form 



5p(r_L,ip,z) 



1 



4(to + 2) 



1 



R ,n+1 



(34) 



with frequencies 



u 



mujj_ 1 + 



k 2 R 2 



(35) 



4m(m + 2) J ' 

Using (M) and (139) in Eq. (E9), we obtain after the evaluation of the expectation values to leading order in (kRj_) 2 



0J+ - OJ- = (oj + - CJ_)o + 



hk 2 



2Mm(m + 2) 2 
-2(m+l) 2 (m + 2) 



m 2 + m + 2 + Am ( — 



2m+2 



)2rn / \ 2m— 2 

+ (m + l)(2m 2 +3m + 2) ( ^- j 



(36) 



where (lo+ — ujJ)q is the result following from (|29j) to which (|36j) reduces for k = 0. For the special case of a centered 
vortex ro = and taking m = 1 the result (36) reduces to a result derived previously in 0. 

Let us now compare these results with those following from an application of the sum-rule approach, developed for 
this problem in Q. As is explained there, the sum-rule approach is based on the assumption that the expectation- 
values of the commutators of certain adequately chosen excitation operators, taken in the ground-state containing 
the vortex, will be exhausted by the two modes ui+, w_ whose splitting is to be calculated. As has been shown 
in H, this approach works very well (i.e., its basic assumption is satisfied to the required accuracy) for the axially 
symmetric case of a centered vortex, where in the Thomas-Fermi and long-wavelength limit the result can be checked 
by perturbation-theory However, we shall now see that for the case of off-centered vortices the sum- rule approach 
fails, because the result of first-order perturbation-theory is not reproduced. 

First, we briefly recall the sum rule approach. Let F± = Yl^-i f±{?j) be the mutually adjoint operators, carrying 
opposite angular momentum, which excite the collective modes ±|m|, respectively. Let be the p-energy weighted 
moments of the dynamic structure factor associated to the excitation operators F± , which can be written as expectation 
values of commutators between the many-body Hamiltonian of the system (H) and F±: 



=<0\[F_,[H,F + }}\0> = 



Nh 2 



<0||V/+| 2 |0> 



(37) 



=<0| [[F^H], {H,F+}} |0>= N <0| j + ] |0> 



(38) 



where |0> may be a vortex-free or a vortex state, and j± — (?i/to)V/± - p. Assuming that the moments are exhausted 
by the modes ±|m|, the shift of the collective frequencies in the TF limit can be calculated as 



h(ui+ — UJ-) = m 2 /rrii 



(39) 



Applying this prescription let us start with the transverse quadrupole modes m — ±2. The excitation operators 
carrying angular momentum m = ±2, are 



f ± = (x± iyf 



(40) 



Evaluating the moments Eqs. (|37]) and (|38|), and from ( |39| ) the following result for the frequency shift for the m = ±2 
modes is obtained J5J: 



2 < l z > 
M <r 2 L > 



(41) 
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Neglecting the microscopic details of the vortex core structure in the density profile, the square radii of the condensate 
can be evaluated using the TF approximation Jll]]. This prescription can be applied to a centered- or non-centered 
vortex. What changes in each of these configurations is the averaged value of the axial component of the angular 
momentum < l z >. Using the results obtained in the previous Section, Eq. ( ^3|) the sum-rule result for the frequency- 
splitting due to an off-center vortex at a distance ro from the symmetry axis, is 



6 



UJ + 



■ Ul- 



Rl 



(42) 



where we recall that the square radius in the transverse direction of a cylindrical condensate in the TF limit is: 



< r, >= 



R 2 4 , 



) 1/2 4 



(43) 



The splitting fl42|) differs from the perturbative result (|33| ) by the factor [1 — (ro/Rx) ] and therefore systematically 
underestimates the frequency-shift, except for a centered vortex tq = 0. 

Let us try the sum-rule approach also on the calculation of the splitting of axial helical vortex modes |^2| with 
low multi-polarity and small wave- number k. The excitation operators of the axial helical vortex modes with wave 
number k and angular distortion m — ±1 and m = ±2 are 



f ± = (x± iy) e 



r_l: / k Z 



(44) 



f± = {x± iy) e 



2 „±ikz 



(45) 



respectively. The calculation of the corresponding moments leads to the following frequency shifts. For the m = ±1 
helical vortex mode: 



k 2 

u + -^ = -<l z > 



(46) 



and for m = ±2 



2 < l z > 
M<r 2 > 



l + ^<ri> 



(47) 



where we have used < r\ >= 3/2 < r\ > 2 . Both expressions agree with the perturbative result for centered vortices 
ro = 0, but fail again to reproduce the correct first-order results for the off-center vortices. 



VI. SUMMARY 



We have calculated the angular momentum of a large cylindrical condensate in the presence of an off-axis vortex 
line and also vortex arrays, within the Thomas Fermi approximation. In addition, the interaction energy of vortex 
lines in the spatially inhomogeneous condensate has been derived for the same case, including the effects of the image 
vortices created by the boundary condition of vanishing normal velocity. It is interesting that the image vortices do 
not contribute to the angular momentum even though they strongly modify the velocity-field around any off-center 
vortex-line. The kinetic energy, on the other hand, is influenced by the presence of the image- vortices. 

We have seen that the contribution of a given vortex to the angular momentum per particle decreases as it moves 
away from the center of the system, which leads to a smooth increase of the angular momentum as a new vortex appears 
at a large distance from the axis of rotation. We have obtained an estimate of the angular momentum corresponding 
to different experimentally observed vortex configurations that qualitatively agree with the continuous behavior 
measured from experimental data || . Finally, we have studied the frequency splitting of low lying collective modes due 
to the presence of an off-centered vortex, including axial vortex-waves, using the perturbative approach jl6| and also 
the sum rule approach A comparison of the results shows, that the simple extension of the sum- rule approach to 
off-axis vortices does actually not work, i.e., the sum-rules are in this case not exhausted by the two nearly degenerate 
modes of interest. 
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Note added: After completion of this work, we have received a preprint p3| reporting numerical calculations of 
critical frequencies and angular momentum of vortex configurations with different vorticity in an elongated trap. The 
vortices numerically obtained in p3|] are not rectilinear but with deformed shape and the longitudinal deformation of 
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FIG. 1. Dimensionless energy of symmetrical vortex-arrays of 1, 2, 3, 4 singly quantized vortices (from bottom curve to top 
curve, respectively) as a function of the rescaled distance of the vortex lines from the z-axis. The condensate has N = 2.5 x 10 
atoms of 87 Rb with R± — 3.8/im, which correspond to the experimental parameters of Ref. M 



FIG. 2. Dimensionless energy of symmetrical vortex-arrays of 1, 2, 3, 4 singly quantized vortices (from bottom curve to top 
curve, respectively) in a trap rotating with frequency S}/2n — 40 Hz as a function of the rescaled distance of the vortex lines 
from the z-axis. The condensate has N — 2.5 x 10 5 atoms of 87 Rb with J?x = 3.8(im, which correspond to the experimental 
parameters of Ref. Q 



FIG. 3. Average angular momentum per particle of symmetrical arrays of 1, 2, 3, 4 vortices (from bottom curve to top 
curve, respectively) in relative equilibrium in a trap rotating with frequency il/2n as a function of Q. The values of M and R± 
are the same as in Figs. 1,2 
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